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I. INTRODUCTION 

The anapole moment was predicted first by Zeldovich 
PQ in 1958 as a new parity violating (PV) and time re- 
versal (T) conserving moment of an elementary parti- 
cle. It appears in the second-order multipole expansion 
of the magnetic vector-potential simultaneously with the 
P- and T— violating magnetic quadrupole moment [2]- 
The nuclear anapole moment was experimentally discov- 
ered in the 133 Cs atom in 1997 [3]. This measurement 
was performed following a proposal by Flambaum and 
Khriplovich [4] , who have shown that the nuclear anapole 
provides the dominant contribution to the nuclear-spin- 
dependent (NSD) parity violating effect in atoms and 
molecules. It can provide important information about 
hadronic weak coupling, which is currently not so easily 
obtained from first-principles nuclear structure calcula- 
tions (e.g., see Ref. [5] and Review [5]). The term in the 
Hamiltonian operator arising from NSD parity violating 
electron-nucleus interaction is 

H A = K NSD ^^—^p(r), (1) 

where knsd is the dimensionless strength constant, 
G F = 2.22249 x 1CT 14 a.u. is the Fermi constant, a 
is a vector comprised of the conventional Dirac matri- 
ces, I is the nuclear spin, r is the displacement of the 
valence electron from the nucleus, and p(r) is the (nor- 
malized) nuclear density. There are three sources for this 
interaction: the first contribution arises from the elec- 
troweak neutral coupling between electron vector and nu- 
cleoli axial- vector currents (V g Atv) [7J. The second con- 
tribution comes from the nuclear-spin-independent weak 
interaction combined with the hyperfine interaction [8]. 
Finally, the nuclear anapole moment contribution scales 
with the number of nucleons, A, with ha ~ A 2 / 3 , and 
becomes the dominant contribution in spin-dependent 
atomic parity violation effects for sufficiently large nu- 
clear charge Z [H [5] . It requires nuclear spin I ^ and 



in a simple valence model has the following value [5] 

k a = 1.15 x 1(T 3 (t^t) A 2 / 3 p l9l , (2) 

Here, K = (-l) /+ iH(/ + 1/2), I is the orbital angular 
momentum of the external unpaired nucleon i — n,p; 
p p = +2.8, p n = —1.9. Theoretical estimates give 
the strength constant for nucleon-nucleus weak poten- 
tial g p rs +4.5 for a proton and \g n \ ~ 1 for a neu- 
tron |10j . The aim of anapole measurements is to pro- 
vide accurate values for these constants. The nuclear 
anapole moment for 133 Cs (7=7/2), containing a single 
valence proton, was measured from the differences in the 
6Sf=4 to 7Sf=3 and 6Sf=z to 7Sf=a hyperfine transi- 
tions as K4=364(62)xl0" 3 [HUTU]. However, the limit on 
9p (ffp = — 2 ± 3 [5]), obtained from the Tl anapole mea- 
surements [TT], seems to contradict the Cs anapole mea- 
surements (g p = 6 ± 1, see [ID]). This indicates that fur- 
ther refinements in the experimental measurements are 
required to obtain high precision results for nuclear spin- 
dependent parity violation effects in atoms. 

In Refs. [T2TTH] it was shown that the nuclear spin- 
dependent parity violation effects are enhanced by a fac- 
tor of 10 5 in diatomic molecules with 2 £i/2 and 2 Hi/2 
electronic states due to the mixing of close rotational 
states of opposite parity (fi-doublet for 2 n 1 / 2 )- Dc- 
Mille and co-workers suggested therefore to measure 
the anapole moment by using diatomic molecules in a 
Stark interference experiment to determine the mixing 
between opposite-parity rotational/hyperfine levels |15j . 
The molecular route opens up the range of systems to be 
studied and should provide data on anapole moments for 
many heavy nuclei. Another motivation comes from a 
possibility to test the standard model. The anapole mo- 
ment contribution is small in light nuclei with a valence 
neutron (Eq. Q). In this case the electroweak contribu- 
tion may be extracted from the measurements of NSD P V 
effects [15 . This contribution has never been measured. 
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We therefore present Dirac-Hartree-Fock (corrected for 
electron correlation) and 4-component density functional 
theory calculations of the electronic factor Wa for the di- 
atomic group-2 and -12 fluorides and a number of other 
diatomic compounds. 

II. COMPUTATIONAL DETAILS. 

For 2 £i/2 and 2 Hi/ 2 electronic states, the interaction 
([!]) can be replaced by the effective operator, which ap- 
pears in the spin-rotational Hamiltonian [HI 115] . 

n A = knsdvva j , (o) 

where S' is the effective spin (discussed below) and n is 
the unit vector directed along the molecular axis from 
the heavier to the lighter nucleus. 

Calculations of the P-odd interaction constant Wa 
were carried out within an open-shell single deter- 
minant average-of-configuration Dirac-Hartree-Fock ap- 
proach (DHF) [16; and within the relativistic density 
functional theory (DFT) |17j , employing quaternion sym- 
metry PH Dl]. We used the DIRAC10 computational 
package [2U] to perform all the calculations. The elec- 
tronic factor Wa is found from evaluating the matrix 
elements of the ap(r) operator in the molecular spinor 
basis [3T]. 2 £i/ 2 and 2 U 1 / 2 open-shell electronic states 
are two-fold degenerate, corresponding to the two pos- 
sible projections of electronic angular momentum along 
n, i.e. = | ± |). When operating within this de- 
generate space, the operator ^=ap(r) is equivalent to 

H^(n x S') (Eq. ([3])). Time- reversal symmetry ensures 
that only the matrix elements that are off-diagonal in 
are non-vanishing. This symmetry rule is encapsulated 
within the effective operator H A S by the angular factor 
(n x S'). Here the effective spin S' generates rotations in 
the degenerate subspace analogously to usual spin oper- 
ator S in a spin- 1/2 system. In the non-relativistic limit, 
S' — > S for 2 £1/2 states. A finite nucleus, modelled by 
the Gaussian charge distribution was employed [55]. We 
note that a nuclear point charge approximation should be 
avoided in calculations of Wa as the resulting singularity 
in the wave function gives unphysical results. 

In the DFT calculations we used the Coulomb- 
attenuated B3LYP functional, (CAMB3LYP*), the pa- 
rameters of which were adjusted by Thierfelder et al. 
[23] to reproduce the PV energy shifts obtained using 
coupled cluster calculations. The newly adjusted param- 
eters are a — 0.20, j3 — 0.12, and \x = 0.90. In order 
to test the stability of the results with respect to the 
choice of the functional, the calculations were performed 
also using the PBE, LDA, and B3LYP functionals. For 
all the systems the Wa parameters obtained using the 
different functionals were in remarkably good agreement, 
within up to ~ 10% (the only exception being RaF, where 
the LDA results were 15% higher than the CAMB3LYP* 
ones). We thus only present the CAMB3LYP* values, 



which are considered to give the best results for parity 
violating properties [23] . 

For the lighter elements (N, O, F, Mg, and CI), uncon- 
tracted aug-cc-pVTZ basis sets were used [24J, [25] . For 
the rest of the atoms, we employed Faegri's dual family 
basis sets [33] . A good description of the electronic wave 
function in the nuclear region is essential for obtaining 
reliable results for parity violating properties [27j . Thus, 
we augmented the basis sets with high exponent s and p 
functions, which brings about an increase of around 10% 
in the calculated values of Wa- The basis sets were in- 
creased, both in the core and in the valence regions, to 
convergence with respect to the calculated Wa constants. 
The final basis sets can be found in the Appendix. 

Where available, we used experimental geometries. For 
molecules where the bond length R e is not known ex- 
perimentally, we optimized the bond distance instead, 
using relativistic coupled cluster theory with single, 
double, and perturbative triple excitations (CCSD(T)) 
[28] . To reduce the computational effort, the Dirac- 
Coulomb Hamiltonian was replaced by an infinite order 
two-component relativistic Hamiltonian obtained after 
the Barysz-Sadlej-Snijders (BSS) transformation of the 
Dirac Hamiltonian in a finite basis set j29[ [30] . Our cal- 
culated R e are typically within 0.01 A of the experimen- 
tal values, where available. The experimental/calculated 
equilibrium distances can be found in Tables [T] and [H} 

Previous investigations for the BaF molecule have 
shown that the electron correlation contribution to the 
Wa constant is non- negligible, and raises its value by 
^20-50% [3U [33]. In this work we use two separate 
schemes to account for the correlation effects: the den- 
sity functional calculations, and correcting the DHF val- 
ues for the correlation contributions using atomic cacu- 
lations, in a manner outlined below. 

The main contribution to the matrix elements of 
the NSD interaction for the valence molecular electrons 
comes from short distances around the heavy nucleus. 
Thus, these matrix elements are strongly affected by cor- 
relations between the core and the valence electrons. 

The total molecular potential at short distances from 
the heavy nucleus is spherically symmetric to very high 
precision; the core of the heavy atom is practically unaf- 
fected by the presence of the second atom. The molecu- 
lar orbitals of the valence electron can thus be expanded 
in this region, using spherical harmonics centered at the 
heavy nucleus, 

\ipv) = a\s 1/2 ) + b\p 1/2 ) + c\p 3/2 ) + d\d 3/2 ) ... (4) 

Only S1/2 and P1/2 terms of this expansion give signif- 
icant contribution to the matrix elements of the weak 
interaction. These functions can be considered as states 
of an atomic valence electron, and are calculated using 
standard atomic techniques in two different approxima- 
tions: one that includes electron correlation and another 
that does not. The correlation factor K tot is then defined 
as K to t = Kw ■ Kei ■ Ke„- Here, Kw is found as the 
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TABLE I. Internuclear distances R e (A) and the P-odd interaction constants Wa (Hz) obtained on different levels of theory: 
DHF, DFT, DHF and DFT scaled for core polariration contribution, ~DFT-K C p and BEF-K C p, DHF corrected for total 
correlation effect, DHF-iC tot , and the final recommended values, taken as W0t(Final)= (W A (DFT)- K C p + Wa(DHF)-K C p))/2. 
Comparison with earlier results is also shown. 
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Semiempirical 
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MgBr 
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2.360 [31] 
1.696 |3T] 
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67 
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Semiempirical 
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164.3 
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Semiempirical 
Semiempirical 


31 
15 




YbF 


Yb 
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35 
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Pb 


2.02fF| 
2.0781511 
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-1382 


3557 
-1349 
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3884 
-1481 


3023 
-1200 
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-1499 


2700 
-720 


Semiempirical 
Semiempirical 




31 
36 
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381 
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2.25^] 


1681 


1465 


2054 


1790 


1773 


1922 


1300 




39 





a Fully relativistic restricted active space configuration interaction method. 

b Relativistic effective core potential (RECP) combined with SCF and an effective operator to account for core-valence correlations. 
c RECP combined with restricted active space SCF approach and an effective operator to account for core-valence correlations. 
d CCSD(T), present calculations 

c RECP combined with spin-orbit direct configuration interaction 

f Quasirelativistic two-component zero-order regular approximation combined with the SCF approach. 



ratio of the matrix elements 



(nsi/ 2 \H A \npi/ 2 ) 

where the matrix element in the numerator includes elec- 
tron correlation, while the matrix element in the denom- 
inator does not. The magnitude of Kyy is larger than 1, 
as the many body corrections due to correlation with the 
core electrons increase the density of the valence electron 
on the nucleus [40], and thus, increase the Wa constant. 

Correlations also modify the expansion coefficients a, 
b, c... in Eq. Q. An estimate of this effect shows that 
it reduces the overall correlation contribution, and pro- 
vides us with Kei and Ke„ , corresponding to the effect 
of correlations on the El amplitude and on the orbital 
energies, respectively. 

Multiplying the result of molecular DHF calculations 
by K to t allows us to include the effect of the most 
important electron correlation contributions (Table I, 
DHF-A'tot). The Appendix contains the details of the 
calculation and the values of K tot for all the atoms stud- 
ied here. 

Core polarization effects are not accounted for in the 
Kramers restricted DFT calculations. These effects are 
included in K to t] however, we also examine their influ- 
ence separately from the rest of the correlation contribu- 
tions. Thus, we define an additional scaling factor, Kqp, 
which only takes the core polarization into account (see 
Appendix for the details of calculation and the values of 
K C p)- Table I contains the DHF and the DFT values, 
corrected for core polarization contribution (DHF- K cp 
and DFT-K C p). 



As the recommended value for the Wa parameter we 
take an average of Wa(DUF)-K C p and W A (DFT) -K C p, 
since we believe that these two values represent our most 
reliable results. 



III. RESULTS AND DISCUSSION 

All the systems under study have a single valence elec- 
tron; the ground state of the PbF molecule is 2 U 1 / 2 , the 
remaining molecules have a 2 Si/2 ground state. Table 
I contains the results obtained for molecules, for which 
earlier calculations were performed. For comparison, we 
present the uncorrelated DHF results, together with the 
different correlation schemes used: the DFT values, the 
DHF and the DFT values corrected for core polariza- 
tion, the DHF values corrected for the overall correlation 
effects, WA(DHF)-K to t, and the final recommended val- 
ues. To the best of our knowledge, these are the first 
DFT calculations of Wa- The spread of the values ob- 
tained by different methods may be used to estimate the 
accuracy of our calculations, as the uncertainty of the re- 
sult is strongly dominated by the correlation correction 
contribution. The difference between the DFT and DHF 
values gives us an estimate for the correlation contribu- 
tion, which is included into DFT and absent in DHF. 
Another (and probably less accurate) method to estimate 
the correlations is to compare Kcp and K tot . From these 
comparisons we have come to conclusion that the accu- 
racy of Wa produced by the metal anapole moments is 
about 15% for E terms and 20-30% for Hi/ 2 terms (such 
as PbF). The accuracy of Wa produced by Br (or other 
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TABLE II. Internuclear distances, R e (A), relativistic fac- 
tors Rw, and the recommended values of Wa constants (Hz), 
taken as WA(Final)= (Wa(DFT)-K cp + W a (DBF) -Ktat))/ 2 
for group-2 and -12 fluorides, and for mercury halides. 



MgF 

CaF 

SrF 

BaF 

RaF 

ZnF 
CdF 
HgF 

HgCl 
HgBr 
Hgl 



12 
20 
38 
56 



30 
48 
80 

80 
80 
80 



Re (A) 

Group-2 

1.750 [31] 
1.967 [3T| 
2.075 [31] 
2.162J3TJ 

2.244M 
Group- 12 

1.766 SI] 
1.99lT 
2.02E 

Mercury 

2.38/F 
2.46^ 
2.73£^ 



Rw 
fluorides 

1.06 
1.17 
1.65 
2.73 
10.34 
fluorides 

1.38 
2.14 
7.03 
halides 

7.03 
7.03 
7.03 



Wa (Final) (Hz) 

5 
11 
52 
153 
1922 

64 
264 
3693 

3647 
3600 
3356 



CCSD(T), present calculations 



halogen) anapole moment is about 30%. 

Most of the previous investigations of Wa rely on 
semiempirical methods, while for BaF, YbF, RaF, and 
PbF ab initio calculations were also performed [32 [33J 
135] [381, 39 . Our results are in good agreement with the 
recent semiempirical values for most of the systems; in 
case of HgF, PbF, and RaF our final Wa constants are 
higher than the previous values by about 30%. 

For BaF we can compare our results to other ab initio 
calculations. Our DHF value is close to that of Nayak 
and Das [32], obtained by a similar method, and to that 
of Kozlov et al. [33], calculated by the combination of 
relativistic effective core potential (RECP) and the SCF 
approach, and corrected for valence-core correlation by 
an effective operator. The value of Wa corrected for cor- 
relation is in very good agreement with both the rela- 
tivistic restricted active space configuration interaction 
(4c-RASCI) result of Ref. [32], and the RECP restricted 
active space SCF approach of [53] , 

Two other systems with ab initio results are YbF 
|35j . calculated using a combination of RECP and 
SCF/RASCF methods, and RaF [35], treated via the 
quasirelativistic two-component zero-order regular ap- 
proximation (ZORA) combined with an SCF approach. 
As both these investigations do not treat electron corre- 
lation, we can compare them directly to our DHF results, 
which are in good agreement. 

Our correlated result for PbF is rather higher than that 
obtained in the previous investigation |38| performed us- 
ing a combination of RECP and spin-orbit direct con- 
figuration interaction, which might be caused by the 4- 
component treatment of relativity in our case. This sys- 
tem is different from the other molecules discussed here, 
due to its ground state. For the 2 Hi/2 electronic state Wa 
vanishes in the non-relativistic limit, since in this limit 
it does not contain the s-wave electronic orbital and can 
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FIG. 1. Scaling of log (j^f) with log(Z) for group-2 and -12 
fluorides 



not provide the matrix element < si/ 2 \Wa\pi/2 >■ The 
effect appears due to the mixing of 2 £i/2 and 2 Hi/2 elec- 
tronic states by the spin-orbit interaction. This gives an 
extra factor of Z 2 a 2 in the Z-dependence of the matrix 
element of Wa in a 2 Hi/2 electronic state. This does not 
lead to a significant suppression in heavy molecules, such 
as PbF; however, in light molecules the matrix element 
of Wa in the 2 n 1 / 2 electronic state is much smaller than 
that in 2 £i/2 state (note that a similar factor Z 2 a 2 also 
makes the interval between the opposite parity A-doublet 
states small; therefore, there is actually no suppression 
in the mixing of these states by Wa)- 

Table [n] contains the recommended Wa constants for 
the group-2 and group-12 fluorides. The magnitude of 
Wa in 2 S 1 / 2 electronic state is expected to scale as Z 2 Rw 



[S], where Rw is the relativistic parameter, 



ft 



2 7 + l 



a B 



3 V2^r A 1 /3 
7 =[l-(Za) 2 ] 1 /2. 



2-27 



[r(2 7 + i)r 



(6) 



-15 



In Eq. ((6), as is the Bohr radius, 7'o = 1.2 x 10 " J m, 
and a is the fine-structure constant. The Rw parame- 
ters are^hown in Table ITT1 for each of the metal atoms. 
In Fig. 1 we plot log ^jj^^ as a function of log(Z) for 

both groups of dimers. For group-2 fluorides the scaling 
is, indeed, close to Z 2 ; however, the interesting feature 
of the plot is group-12 fluorides, where the Z-dependcncc 
is much more advantageous, of Z 2 5 . This is due to the 
filling of the lower lying ci-shell, which expands relativis- 
tically and thus increases the effective nuclear charge, 
leading to an enhancement of relativistic and PV effects 
[42], and an increase of Wa- It should be noted that 
measurements for Wa in MgF and CaF may be used to 
test the standard model, since the anapole moments of 
Mg and Ca are small and the electroweak contribution 
to NSD PV electron- nucleus interaction is important. 

Table [TT] also shows the Wa constants of mercury 
halides. The Wa values are very close for all the halo- 
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gens; the situation is similar for zinc and cadmium 
halides. However, molecules containing heavier ligands 
might have an experimental advantage. Due to the higher 
reduced mass, the interval between the opposite parity 
rotational levels becomes smaller, and thus larger PV ef- 
fects can be expected, and a smaller magnetic field would 
be needed to reduce the interval between the levels. 
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IV. APPENDIX 

A. Basis sets 

Table Mil contains the basis sets used in the calcula- 
tions of the Wa constants. For the lighter elements (N, 
O, F, Mg, and CI), uncontracted augmented correlation- 
consistent valence triple-^ (AVTZ) Gaussian basis sets 
were used [Ml Hi]. For the rest of the atoms, we employed 
the Faegri's dual family basis sets [53], augmented with 
higher orbital momentum and diffuse functions. Good 
description of the electronic wave function in the nuclear 
region is essential for obtaining reliable results for parity 
violating properties [57] ■ Thus, we also augmented the 
basis sets with high exponent s and p functions, which 
brings about an increase of around 10% in the calculated 
values of Wa- The basis sets were increased, both in the 
core and in the valence regions, to achieve convergence 
with respect to the obtained Wa constants. The same 
basis sets were employed in the relativistic CCSD(T) ge- 
ometry optimizations (where experimental geometry was 
unavailable). However, here we leave out the high ex- 
ponent s and p functions, as these contribute little to 
molecular geometries. 



B. Calculation of the correlation factors Kcp and 

K tot 

The scaling factors Kc p and K to t are used for estimat- 
ing the effect of the core polarization and of the overall 
correlation effects on Wa for the molecules. These fac- 
tors are found using two sets of atomic calculations, one 
that neglects electron correlation and on that includes 
it. The calculations, and the derivation of the Kcp and 
K tot factors are described below. 



1. Atomic calculations without electron correlation 

We use the relativistic Dirac-Hartree-Fock method 
(DHF) to perform atomic calculations [UJ. I n atomic 



TABLE III. Basis sets employed in the calculation of the Wa 
constants. All elements with Z > 17 are described by the 
Faegri basis sets [26] augmented by high exponent, diffuse, 
and high angular momentum functions. 



Atom 


Z 


Basis Set 


N 


7 


aug-cc-PVTZ 





8 


aug-cc-PVTZ 


F 


9 


aug-cc-PVTZ 


Mg 


12 


aug-cc-PVTtf] 


CI 


17 


aug-cc-PVTZ 


Ca 


20 


20sl8p9d6/l# 


Zn 


30 


21sl9pl0d7/2g 


Br 


35 


21s20pl0dl0/lp 


Sr 


38 


21s20pl2rf9/2g 


Zr 


10 


21s20pl2rf9/2g 


Cd 


18 


22s20pl2rf9/2g 


I 


53 


22s21pl2rfll/2t/ 


Ba 


56 


24s22pl5dl0/2s 


La 


57 


24s22pl4dl0/2ff 


Yb 


70 


26s21pl4dl0/23 


Hg 


80 


25s21pl5<210/2g 


Pb 


82 


25s22pl6dl0/2s 


Ra 


88 


26s23pl6dll/2# 



augmented by 2 high exponent s and 4 high exponent p 
functions. 



units (|e| = 1,H = l,m e = 1), the single electron DHF 
Hamiltonian is given by 

# = ca-p + (/3-l) C 2 - - + V e (r), (7) 

r 

where a and (3 are the Dirac matrices and V e (r) is the 
self-consistent DHF potential due to atomic electrons. 
In order to take into account the specifics of diatomic 
molecules, we use a slightly modified V e (r) potential com- 
pared to standard atomic techniques: 



Ve(r) = V, 



DHF 



1 



(8) 



Here Vdjj^ " is the self-consistent DHF potential of the 
closed-shell core of the heavy atom, N is total number of 
electrons in this atom (N = Z for a neutral system), N v 
is the number of valence electrons, r is the distance to the 
heavy nucleus, and R e is the distance between the nuclei 
in the molecule. The second term in ^ represents the 
spherically symmetric contribution from the valence elec- 
trons that are assumed to be moved to the second atom. 
Its form is chosen to have the correct — 1/r asymptote 
for the total potential at large distances in the case of 
neutral molecule. 

The self-consistent DHF procedure is first done for the 
ion, from which valence electrons are removed. Then the 
core potential V D - S ^ r " is frozen and valence S\/2 and P1/2 
states are calculated by solving the DHF equations for 
the valence electron 



(H - e v )i>v = 0, 



(9) 
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where H is given by (|7| and (|8j). 

For example, in the case of the BaF molecule, Z — N — 
56, N v = 2, and R e = 2.162 A [3"T]. The self consistent 
DHF procedure is performed for the Ba 2+ ion. The Qsi/ 2 
and 6pi/2 valence states are calculated in the potential 



56 

V(r) = +y D H F (Ba 

r 



2+\ 



1 



y/r 2 +R 2 e 



(10) 



These 6s 1/2 and 6pi/2 states are used in the denominator 
of(§. 



U. Inclusion of electron correlation 

We include two important classes of electron correla- 
tion corrections: the core polarization correction and the 
Brueckner-type correlations. These types of correlations 
dominate in S1/2 and P1/2 atomic states and their inclu- 
sion leads to an accuracy of a few percent for the matrix 
elements [43] . 

The core polarization can be understood as the change 
of the self-consistent DHF potential due to the effect of 
the extra term (the weak interaction operator Ha) in 
the Hamiltonian. The inclusion of the core polarization 
in a self-consistent way is equivalent to the well-known 
random-phase approximation (RPA) (see, e.g. [13]). The 
change of the DHF potential is found by solving the 
RPA-type equations self-consistently for all states in the 
atomic core. The RPA equations have a form of the DHF 
equations with the right-hand side: 



(H - e c Wc = {-H A + SV A )^ C 



(11) 



where Hq is the DHF Hamiltonian (|7|), index c enumer- 
ates states in the core, Sipc is the correction to the core 
state c due to weak interaction Ha, and 5Va is the cor- 
rection to the self-consistent core potential due to the 
change of all core functions. Once SV A is found, the core 
polarization can be included into a matrix element for 
valence states v and w via the redefinition of the weak 
interaction Hamiltonian, 



(v\H A \w) ->■ (v\H A + SV A \w) 



(12) 



We then obtain the scaling parameter for core polar- 
ization, Kcp, from 



K, 



CP 



<e^Ai^r /2 > 



(13) 



The values of Kqp for all the atoms under study are 
presented in Table [TV] 

In contrast to the core polarization correction, which 
can be reduced to the redefinition of the interaction 
Hamiltonian, the Brueckner-type correlations can be re- 
duced to the redefinition of the single-electron orbitals, 
by replacing the DHF orbitals by the Brueckner orbitals 



(BO). Brueckner correlations describe the interaction be- 
tween the valence and the core electrons. These correla- 
tions can be included with the use of the so-called corre- 
lation potential £, which is defined in such a way that the 
average value of £ over a valence state v is the correlation 
correction to the energy of this state: 



Se„ = («|£|u) 



(14) 



The correlation potential £ is a non-local operator sim- 
ilar to the DHF exchange potential. It can be calcu- 
lated by means of the many-body perturbation theory 
(MBPT) in the residual electron-electron Coulomb inter- 
action. The expansion starts from second order in this 
interaction and in most cases this is the leading term. 
We use the B-splines in a box [33] and the second-order 
MBPT to calculate £. The Brueckner orbitals for the va- 
lence states are found by solving the DHF-like equations 
with an extra operator £ included: 



(H + £ - e v )^° = 



(15) 



Solving these equations gives new energies and new wave 
functions for the valence states. For all atoms considered 
in present paper the inclusion of the second-order cor- 
relation potential £ leads to a few percent accuracy for 
the energies of the Si/ 2 and P1/2 states. Matrix element 

of the operator Ha between valence states v and w, in 
which both type of correlations are included, is given by 



The K w factor (Eq. ^) is then reduced to 



K 



w 



(^° /2 \H a + SVa\^P Pi/ 



(16) 



(17) 



An additional effect of correlation, not taken into ac- 
count in the above is the change in the expansion coeffi- 
cients a, b, c... in Eq. Q. In order to treat this effect 
we turn to the following expression, 



(nsi/2 \H A I npi /2 ) (npi/2 1 El | ns 1/2 ) 



E,, 



'1/2 



E 



(18) 



"Pl/2 



Such an expression appears in atomic parity violating 
elecromagnetic amplitudes and, in addition to the weak 
matrix element, it contains the El amplitude and the en- 
ergy denominator between the s and the pi/2 states. In 
molecules, use of this expression may be justified by the 
ionic bond model, where the electron that moves to the 
halogen polarizes the metal atom and produces a mix- 
ture of s and p orbitals. Table [TV] contains the obtained 
K\y (corresponding to (ns\H A\ n p)) , along with Kei and 
K.E n factors, which describe the effects of core polariza- 
tion (Eq. plj)-([l2"|)), and the correlations (Eqs. pl|- 



(IT)) on the El amplitude and the energy denominator, 
respectively. The final rescaling parameter K tot is the 
product of all three factors 



Ki, 



Kw ■ Kei ■ Ke„ 



(19) 
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TABLE IV. KCP, K w , K E \ , K Erl , and K tat correlation scal- 
ing vactors for the atoms under study. 



Atom Kqp 


K w 


K E i 


K En 


Ktot 


Ms 


1.24 


1.38 


0.97 


0.96 


1.29 


Ca 


1.28 


1.54 


0.91 


0.90 


1.26 


Zu 


1.12 


1.38 


0.87 


0.91 


1.09 


Br 


1.31 


2.41 


0.84 


0.80 


1.34 


Sr 


1.26 


1.60 


0.88 


0.86 


1.21 


Zr 


1.20 


1.40 


0.88 


0.91 


1.12 


Cd 


1.11 


1.40 


0.80 


0.86 


0.96 


Ba 


1.26 


1.69 


0.84 


0.82 


1.16 


La 


1.22 


1.54 


0.85 


0.86 


1.13 


Yb 


1.20 


1.63 


0.83 


0.83 


1.12 


Hg 


1.09 


1.33 


0.74 


0.86 


0.85 


Pb 


1.10 


1.22 


0.77 


0.95 


0.89 


Ra 


1.22 


1.64 


0.91 


0.81 


1.21 



The final K to t factors can be found in Table |IV| and 
are used to scale the DHF Wa parameters to include the 
effects of electron correlation. (Note that that for PbF 
we have calculated only correlations between the valence 
electrons and the Pb core; the correlations between the 
valence electrons 6s and 6p are not included and should 
be treated separately using a different technique). The 
remaining correlation corrections, which are often called 
the weak correlation potential 5Sa or structural radia- 
tion, are suppressed by a small ratio e v /e c ~ 0.1, where 
e v is the valence electron energy and e c is the core elec- 
tron energy. The effect of 5Sa usually does not exceed a 
few percent [13 SS] . 
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